Statistical analysis of complex systems with nonclassical invariant measures 
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I investigate the problem of finding a statistical description of a complex many-body system whose 
invariant measure cannot be constructed stemming from classical thermodynamics ensembles. By 
taking solitons as a reference system and by employing a general formalism based on the Ablowitz- 
Kaup-Newell-Segur scheme, I demonstrate how to build an invariant measure and, within a one 
dimensional phase space, how to develop a suitable thermodynamics. A detailed example is provided 
with a universal model of wave propagation, with reference to a transparent potential sustaining gray 
solitons. The system shows a rich thermodynamic scenario, with a free energy landscape supporting 
phase transitions and controllable emergent properties. I finally discuss the origin such behavior, 
trying to identify common denominators in the area of complex dynamics. 
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I. INTRODUCTION 

The field of Statistical Mechanics (SM) delves with 
the problem of deducing macroscopic properties of mat- 
ter stemming from microscopic states, specified by 
mechanical-like phase spaces whose dimensionality de- 
pends on the atomic degrees of freedom of the system. 
Classical SM lies its foundations in the orthodicity prob- 
lem posed by Boltzmann in the XIX century and in 
particular in the search of invariant measures [2| whose 
infinitesimally change leads, in the thermodynamic limit, 
to thermodynamic relations among internal energy, vol- 
ume, pressure and average kinetic energy. Boltzmann's 
work leads to three models of thermodynamics, namely 
the microcanonical, canonical and grand canonical en- 
sembles, each characterized by a diverse invariant mea- 
sure of the phase space In more recent times, the 
same invariant measures have been found by Shannon 
in a famous paper by employing a rather different 
approach based on an variational minimization problem 
connected to an information functional defined axiomat- 
ically. Boltzmann and Shannon's approach has in com- 
mon the search of invariant measures under specific con- 
straints (i.e., the thermodynamic relations or the mini- 
mization of the information functional), whose role is to 
define models of thermodynamic with the minimum in- 
formation possible on the microscopic state. However, a 
natural question therefore arise, and concerns the possi- 
bility to develop a SM description of microscopic mechan- 
ical systems whose invariant measure is not constrained 
by any a-priori condition. This problem acquires funda- 
mental importance in the field of Complex Systems (CS). 
These are characterized by the presence of a large number 



* landrea.fratalocchi@uniromal.itl | www. primalight . org| 



of interconnected or interwoven parts, which act together 
to sustain emergent properties of the whole ensemble Q . 
Illustrative examples of CS are found in many different 
fields including — and not limited to — neural networks, 
colloids, complex liquids and glasses 0-0] • Contrary 
to classical thermodynamic systems, mainly described 
within two scales (micro and macro), complex systems 
exhibit a full set of intermediate length scales each char- 
acterized by specific properties (see e.g. [T^ in the con- 
text of structural glasses). In this article I focus my anal- 
ysis on the so called mesoscopic scale Q , which contains a 
large (but finite) number of atomic parts n. This length 
scale allows for a general treatment of the problem as 
it does not require any hypothesis on the structure of 
the invariant measure at the thermodynamic limit, when 
n — >■ oo. With these premises, my initial question be- 
comes: with reference to a dynamical many-body ensem- 
ble, whose invariant measure is not constrained by any a- 
priori condition, could a statistical mechanic framework 
be developed at the mesoscopic scale? 
An important benchmark for such problem is provided by 
the theory of solitons, which are exact, localized particle- 
like solution of nonlinear systems integrable through the 
inverse scattering transform jlh . |12| . Although there are 
well documented observations of soliton waves in the past 
centuries [l^ , their mathematical description has become 
clear only in recent times Owing to their robust 

particle nature, solitons are preferential carrier of energy 
transport and are therefore ubiquitous in physics |15l424l | . 
According to our initial question, we ask what is the in- 
variant measure of a mesoscopic many-body ensemble of 
solitons, its relationship with canonical ensembles and, 
more important, how to develop a proper statistical me- 
chanics framework to describe the appearance of (any) 
emergent dynamics. The development of such theory is 
not only of interest to the field of mathematical physics 
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and statistical mechanics, but could be the backbone for 
novel studies aimed at describing extreme events occur- 
ring in nonlinear dynamics, such as e.g., the formation of 
rogue waves, which are attracting a large scientific inter- 
est [Uli^. 

This paper is organized as follows. Section HIl presents a 
general formulation of the problem within the Ablowitz- 
Kaup-Newell-Segur scheme. An invariant measure of the 
soliton system is calculated in Sec. Ill Al while its thermo- 
dynamics analysis is developed in Sec. IIIBI The general 
theoretical framework is then applied to a specific case in 
Sec. mil encompassing both a thermodynamic and time- 
domain analysis. Sec. IIVI finally, discusses the origin 
of the solitons complex behavior, highlighting similari- 
ties and differences with other systems showing complex 
dynamics. 



II. GENERAL FORMULATION THROUGH THE 
AKNS SCHEME 

In order to pursue a general theory, I eroploy the 
Ablowitz-Kaup-Newell-Segur (AKNS) scheme [27| whose 
general framework encompasses a large selection of inte- 
grable equations in both (1 -f- 1) and (2 -I- 1) dimensions, 
including the Korteweg-de Vries (KdV) , the Kadomtsev- 
Petiashvili (KP), the Sine-Gordon (SG) and Nonlinear 
Schrodinger (NLS) equations. I begin by introducing the 
pair of equations: 



dx 
d_ 

dt 



- U 
-A 



(1) 
(2) 



with F = 
matrices: 
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1/^2 



a two-element vector and the 2x2 
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A = 
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being A the spectral parameter, q{x,t), r(x,t) generic 
complex potentials and Aij = Aij(x, t, X) unknown func- 
tions of the spectral parameter A. The compatibility be- 
tween Eqs. ([Ij and ([2|) yields the zero curvature condi- 
tion: 



d_ 

dx 



A- 



d_ 

dt 



U + [A, U] = 0, 



(4) 



which gives an infinite hierarchy of integrable nonlin- 
ear partial differential equations generally expressed as 
f{q,r) = 0. To ensure the existence of F on the infinite 
line X G [— cxD,oo], I need to specify the behavior of the 
potentials q and r at the boundaries |a:| — » 00. Following 
the original paper of Ablowitz et al. [27| . and without 



loss of generality, I assume that both q and r decay suffi- 
ciently rapidly to (we will see in Sec. IIIII how to handle 
the finite density case, i.e., g,r -> pe*"^ for x -> ±00). 



A. Looking for an invariant measure 

The solution F can be expressed by a linear combi- 
nation of the Jost functions F±(x, A), whose asymptotic 
behavior is found by solving ^ for — >■ 00: 



F+{+oo, A) = (Toe 



00, A) = a^e 



(5) 



being ai {i = 1,3) one of the SU{2) generators (ctq is the 
identity matrix): 



CT2 



1 

1 



1 



0-1 



1 

1 

1 ■ 

-1 



(6) 



Each Jost function F± provides a basis for the solution 
of Eqs. ([l}, hence, they are linearly dependent and con- 
nected via the scattering matrix T(A): 



F^=F+T{X), 



T = 



a{X) b*{X) ' 
b{X) -a*(A) 



(7) 



The knowledge of T yields an equivalent representation 
of the dynamics f{q, r) = via the spectral transform S: 



S{q,r) = 



6(A„) 



da(A„) 



(8) 



with —00 < A < 00 and n — 0, N . The transform S is 
composed of: 

i) the reflection coefficient R, which composes the con- 
tinuum spectrum describing radiation; 

ii) the discrete eigenvalues A„ [originated from the zeros 
of a(A)] and the normalization constants c„. This consti- 
tutes the discrete spectrum defining solitons. 

Thanks to the spectral transform, the dynamics of both 
q and r is expressed in terms of the time evolution of the 
spectral coefficients a and b. From the zero curvature 
condition (|4]), straightforward calculations [13, yield: 



a(A, t) = a(A, 0) = a(A), 6(A, t) = 6(A, 0)e" 



(9) 



being a(X) a generic function whose form depends on the 
specific equation considered. According to Eqs. the 
scattering coefficient a(A) is a constant of motion and 
can be therefore employed to derived an invariant mea- 
sure for the overall system. This can be obtained from 
the soliton discrete eigenvalues Ai, which constitute the 
time- invariant part of the spectral transform (|H]). More 
specifically, considering the one dimensional phase space 
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given by the infinite line A e [—00,00], I define the soh- 
ton invariant measure /i(A) from the eigenvalues Density 
Of States (DOS): 

with Xj being a soliton eigenvalue. The advantages of 
are twofold: 

i) it yields an exact description of the dynamics of soli- 
ton waves, as they originate from the discrete spectrum 
which depends on A„; 

ii) it is defined on a single dimensional phase space. This 
allows for a simpler analysis with respect to the time de- 
pendent description of q, r, which gives rise to infinite 
dimensional phase spaces harder to handle. 

It is important to stress that the invariant measure ([10]) 
exists only in nonlinear systems supporting soliton waves. 
In linear regime, in fact, the spectral transform reduces 
to the reflection coefficient R only [29]. This quantity is 
not an invariant of motion (see Eqs. |9]), and no spectral 
invariant measure can be built in this case. 
In summary, a probability function based on the con- 
served part of the spectral transform is an invariant mea- 
sure for an ensemble of solitons, even in the presence of 
radiation, which allows for the development of an exact 
system description within a one dimensional phase space. 



B. Mesoscale thermodynamics through Lie 
transformation groups 

Equations ©-(HU]) possess the following two proper- 
ties: 

i) Solitons eigenvalues A„ are arbitrary in the interval 
A e [— oo,-|-oo]. The inverse scattering, in fact, guaran- 
tees a 1 : 1 relationship between a chosen distribution of 
soliton DOS and a scattering potential V'Oj which may be 
calculated by solving the corresponding Gelfand-Levitan- 
Marchenko integral equations [ij. In consequence of 
that, the soliton invariant measure (|10p does not satisfy 
any entropy maximization principle and could not be ex- 
pressed in terms of the three classical thermodynamic 
ensembles; 

ii) Any invariant of motion, including the Hamiltonian, 
follows the single-soliton factorization: 

'H = "Hn ; 'Hn='Hn{\n), (H) 

n 

which stems from the factorization property of the spec- 
tral transform: 

N 

S = {n = 0, ...,iV, A„,c„} = ^5„, Sn = {A„,c„}, 

(12) 



and yields a microscopic mechanical system (a soliton 
gas) described by a purely kinetic Hamiltonian: 

n n 

It is worth to observe that the wave-particle "dualism" 
originated by Eqs. (|4]) and ([TS]) — i.e., a wave field 
[q{x,t),r{x,t)] associated to particle modes Xj with mo- 
menta pj — goes beyond the ordinary wave-particle du- 
alism of classical optics or quantum mechanics, and is 
only manifested in nonlinear equations exhibiting soliton 
waves. In quantum mechanics, for example, due to the 
Heisenberg uncertainty principle it would never be possi- 
ble to observe a double localization — i.e., the appearance 
of position (momentum) eigenmodes which are localized 
in momentum (position) — or localized wavepackets able 
to overwhelm spreading effects. For soliton waves, con- 
versely, we do observe an ensemble of kinetic eigenmodes 
with specific momenta p„ = y/2'Hn, whose spatial form 
remains localized in space. This is a peculiar property of 
solitons and stems from their pure nonlinear nature. 
In order to provide a statistical mechanics analysis of a 
solitons ensemble, we come back to our original question 
of how to describe of a microscopic mechanical system un- 
der a non canonical measure. I address this problen i by 
mutuating ideas from the thermodynamics of chaos [30] ■ 
In classical statistical mechanics a special role is deserved 
to canonical probability distributions, as they form the 
only ensemble (canonical) that satisfies the thermody- 
namics relations even outside the thermodynamic limit. 
Another important property of canonical distributions, 
which I will exploit here, is their algebraic structure. I 
begin by considering the soliton DOS ([TO)) . defined on 
the line interval A G [Xmin, Xmax], and partition the lat- 
ter into — {Xmax — ^min)/^ boxes of size e. I then 
define pi as the probability of finding an eigenvalue in 
the z— th box: 

P^ = / /i(A)dA, 1 = 1,2,..., TV,. (14) 

"'A„i„+(i~l)e 

Then, the application of a canonical measure on pi [a 
special case of Escort Distributions (ED) [30j]: 

P^^f{p^,P)^^^, (15) 

yields a Lie group [sil, [sl] of transformations acting 
on the space Pi, with the infinite dimensional param- 
eter of transformation /3 = l/T playing the role of 
an inverse temperature. The identity element is ob- 
tained for f3 — I, while straightforward calculations yield 
the Lie group composition law function 0, which reads 
^(/3i)/32) = /3i/32- Such a Lie group of transformations 
can be linked to the field of statistical mechanics by defin- 
ing the following partition function: 

Z(/?) = ^pf =exp(-*). (16) 

j 
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Thanks to this definition, in fact, the representative prob- 
abihty Pi becomes: 



P, = exp(^' - pS,), 



(17) 



with Si = —In Pi acting as the energy of the i-th mi- 
crostate. The Helmholtz free energy J- is then given by 
the corresponding expression in the canonical ensemble: 



(18) 



To develop a mesoscale thermodynamics, I study how 
the Free-energy approaches the thermodynamic limit, de- 
fined by 



lim 



V 



with V — — Ine being the volume. For e 
sponding partition function scales as: 



(19) 

0, the corre- 



(20) 



It is worth to observe that the calculation of the rescaled 
Free Energy x does not need the real computation of the 
thermodynamic limit ^ — > oo, but is obtained from the 
scaling of the partition function as the volume increases. 
In conclusion, in the case of a mechanical phase space 
possessing an arbitrary measure, a statistical description 
is possible by exploiting the Lie group structure of the 
canonical ensemble. This yields a general mesoscale ther- 
modynamics framework [i.e., Eqs. ([T8|) - ([20l) ]. In order to 
illustrate these concepts, I will consider in the next sec- 
tion a specific example. 



III. APPLICATION TO THE NONLINEAR 
SCHRODINGER EQUATION 

Among the various integrable equations described by 
the AKNS scheme, I discuss a universal model of non- 
linear wave propagation in dispersive media, namely the 
NLS equation: 

on the infinite line x G [— oo, oo] with the general bound- 
ary condition: 

4> pe*"^, X ±oo. (22) 

The NLS arises from the AKNS system by taking: 



1 



dx dx 



l(7i 



0-2 



\ dx 



dip 
dx 



XU. 



(23) 



The universality of this model [l5| makes it extremely 
interesting on the experimental side [l^ [33l - l40| , while its 
specific nonlinear sign (defocusing) allows to extend the 
application domain of this thermodynamics framework 
to the general boundary conditions expressed by (|22)) . 



A. Spectral transform of a reflectionless scattering 
potential 

A refiectionless potential is, by definition, an input po- 
tential = "Pi^i 0) whose spectrum contains just the 
discrete part (i.e., soliton waves). For the defocusing 
NLS, I consider the input: 

tpi^i = — p{w + iv)[w ia,n\\{wx) + iv]^ w'^+v^ = \. (24) 

When p = 1, Eq. (|24)) yields the one-soliton solution of 
the NLS with "grayness" w and velocity 2v [4l|, ll^j- In 
the following, I will consider the general case p ^ 1 and 
find the conditions to have a refiectionless potential. I 
begin my analysis with the spectral problem: 



dF 

— = [7(^-0, x,A)F, 



(25) 



also known as linear auxiliary problem [4l|, and in par- 
ticular by finding its Jost solutions F±. These are the 
solutions of the following system: 



^ofTi + ^^asVi, Voe3?, (26) 



dF± 
dx 



being F±(x, A) a 2 x 2 matrix with asymptotic conditions 
F±(±oo,A) ~ F±ao, the latter calculated from ([26l) at 
X — ±oo: 



F^r 



1 



1 



+00 ^ ^ —OO: 



(27) 



being k = -^/A^ — 4p2 



Case V — 



When V = 0, the input ipQ turns into the so called 
tanh-like potential, which is known to be integrable [42| - 
131 ■ Nevertheless, a complete derivation of the spectral 
transform is missing in the literature. For this reason, I 
will provide a thorough analysis of this limiting condition 
in the following paragraphs. This derivation will be also 
helpful when considering the more general case v ^ 0. 
I begin by calculating the Jost solution F_. I apply the 
transformation: 



F_(a;,A) = F_ooS(x,A), 



(28) 



5 



with E being a 2 x 2 matrix of components: 
E 

and asymptotic value 



Ell Ei2 
E21 E22 



E{-oo,X) 



1 
1 



(29) 



(30) 



thus obtaining a hnear auxiliary problem for E: 

= V^[2W3 + A.,e— (31) 

By employing the change of coordinate 

1 + tanh X 
s = 2 ' ^^^^ 

I transform (j3ip into an Hypergeometric equation for 
Eii: 



s(l - s) 



s + ^^)^ + p'Ell=0, (33) 



defined by the following Riemann V symbol: 

r 1 00' 
Eii^Vl -p 

U+^l -*f P 

The symmetry properties of pi^ : 

E21 = Ell 



k—^ — k 7 



and: 



E22 — Ell, 



E12 — E21 , 



(34) 



(35) 



(36) 



allows to express all the elements of the unknown matrix 
E through the single term En defined by ([M)) . In par- 
ticular, from the properties of Hypergeometric functions 
near s — [4^] (i.e., for x —00): 



E = 



2Fi[ - p,p,i^,s 



l-/9 + z|,l + p + zf,2 



2 2 , S 



2 • C 



FAl- p-i^,l + p 



i- 2 - 



i 2 , s 



■2 • c 



2F1 



(37) 



As the reader can verify by straightforward algebra, in 
the limit s — 0, E aQ. The constant c is to be deter- 
mined from the prolongation of F^ to a; — >■ -|-cxi, needed 
to relate F^ to As stated in ([7]), the Jost solutions 
F± are not independent and are related through the lin- 
ear scattering matrix T(A). The latter, considering our 
boundary conditions, is expressed as [41| : 



F_ix,X)^F+{x,X)Tp{X), 
a{X) b*{Xy 



TpiX) 



b{X) a*{X) 



(38) 



with the transition coefficients a(A) and b{X) related by 
|ap — |6p = 1. From the latter condition and through the 
analytical continuation of the Hypergeometric functions 
near s = 1 (x — )■ -l-oo), I obtain the the constant c = 
2iXp/k{2i + k) and the transition coefficients: 



6 = 



'A.r(^ + p)r(i--p)' 



2r(i-p)r(p)- 



(39) 



The coefficients a and b of the scattering matrix allow to 
calculate the spectral transform iSjV'o} of the real input 



pulse ipQ. In particular, for integer values of p: 
R{X) = 0, 



A±„ = ±2Vp2- (p-n)2, n = 0, . . . , M (40) 
while c„ can be calculated from [4l|: 

Cn = ^. (41) 

a(x,, + z^4p^-Xl] 
and reads (see App. |A|) : 
c±„ = (-l)"+i=^M!i!LZfl) , n = 0,...,M (42) 

llj/±r, (s±n - Sj) 

with s±n = X±n + '2i{p — n), j £ [~M,AI] and M — 1. 
The spectral transform Slipo}, for integer p, is therefore 
expressed as: 

ri?(A) = o, 

S{^o} = {^±n = ±2v/p2 - (p-n)2, (43) 

- c pn+i n,(«±^-«;) 

- ^ ll,v±-^(«±"-''^-)' 
with n = 0, . . . , p — 1. 
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FIG. 1. (Color Online), (a) Theoretical (solid line) and nu- 
merical (red circles) soliton invariant measure versus A/2p 
calculated for grayness w = 0.9; (b) Theoretical (solid line) 
and numerical (red circles) free energy X versus inverse tem- 
perature P for w = 0.9. 




2. General case v ^ 

By applying the coordinate change s — (1 + 
tanh?i;a;)/2, with the F-Homotopic transformation 
= (s - \)-^^s-^i^^, I reduce Eq. ^ to the 
normal form of Heun equation (46j : 



1 c oo 
ip = V { Q a s q 

1-7 2 P 



(44) 



FIG. 2. (Color Online). Free energy profile for varying tem- 
perature and grayness w. 



B. Mesoscale thermodynamics 

The soliton invariant measure, for the scattering po- 
tential '00 J reads: 



p-i 

E 

rn—l — p 



(46) 



with c = {w—iv) 121X1, j — d — l-~ik/2'w, /3 = —ik/2'w+p, 
a — —ik/2w — p, and q being the accessory parameter 
whose cumbersome expression will not be reported here. 
The Jost solutions and the a and b coefficient of the scat- 
tering matrix can be then constructed by applying the 
same analysis of Sec. IIII A 11 However, the connection 
problem associated with the Heun equation, i.e. the ana- 
lytic prolongation of Heun functions, is to date unsolved 
[46| . A complete derivation on the spectral transform is 
then not possible in this case, due to the impossibility of 
calculating the zeros of a and b in analytic form. How- 
ever, closed form expression for the radiation spectrum 
and the soliton eigenvalues are accessible by exploiting an 
original expansion with Hypergeometric functions [47|. 
In particular, for integer p: 



having ordered the eigenvalues (|45p in ascending order. 
At the mcsoscopic scale p ^ 1, I further simplify the 
soliton DOS by introducing the variable y — {\ — n/ p)w, 
continuous in the range [0, w], and then rescale the soliton 
eigenvalues A — > A-2p, with |A| € [-\/l — w"^, 1]. With this 
position, for A > 0, Eq. becomes: 



p(\) 



A 



2wVl - A2 



(47) 



This analytical expression compares well with numerical 
estimates (Fig. [T]-a). As seen from Eq. (|47l) . the soliton 
invariant measure depends on two parameters: the soli- 
ton eigenvalue A and the grayness w; the latter condition 
adds one more coordinate to the free energy x = x(/^i 
which becomes function of two different variables. 



fi?(A)-0, 

S{M = < A±„ = ±2^p^ -[p-nYw\ (45) 
[Aq = 2vp, 

with n=l,...,/3 — 1. 

The scattering potential V'Oj expressed by Eq. (p4|) . is 
therefore a reflectionless potential for integer amplitudes 
p, and contains 2p — 1 soliton particles whose eigenvalues 
are expressed by 



1. Free energy x versus temperature: Emergence of pfiase 
transitions 

I begin by studying the structure of the free energy 
landscape for a varying temperature I//3, thus analyzing 
the phase transition scenario of the system. In order to 
apply the thermodynamic framework discussed in Sec. 
IIIBl I first partition the spectral gap [—1,1] into — 
2/e boxes of length e and then calculate the probability pi 
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of finding an eigenvalue in the box [— 1 + (i — l)e, —1 + ie]: 

p-i 



Pi = / ^(A)dA 



2p 



2p 



n— 1— p 

dA(5(A — A„) oc 



(48) 



where in the last step Eq. (|T7)) has been used. When 
e — >■ 0, two principal scaling regimes in (|48p emerge. In 
particular, in the case of |A| ~ 1, an expansion of pi for 
I « yields at leading order a factor proportional to -y/e, 
while for |A| ~ 0.5 the leading order of the expansion 
becomes of 0(e). In summary: 

{Pj oc e for |A| ~ 0.5, 
Pj oc for |A| fa 1. 

The partition function Z then scales as: 

Z = J2Pj -iNe-2)el^+2e^/^, (49) 

with two boxes near A = ±1 scaling like e^/^, and the 
remaining — 2 scaling like . In the limit e — > 0, the 
partition function reads: 



(50) 



This implies X — mm[/3 — 1, /3/2] and leads to: 

x = i':, ^' '1-':^ (51) 



I, /? > /3c 



being /3c = 2 a critical point given by /3c — 1 = /3c/2. The 
free energy (|5ip . being continuous but not differentiable 
at (3c, supports a first order phase transition as the in- 
verse temperature (3 decreases below Pc — 2. As seen in 
Fig. (HJb), this result is in agreement with the numeri- 
cally calculated free energy X from the soliton invariant 
measure 



2. The complete scenario: controllable emergent properties 
and free-energy collapse. 

The effects of the grayness parameter w on the dy- 
namics are then evaluated by inspecting the soliton phase 
space structure. In particular, a decreasing oiw from 1 to 
is accompanied by the shrinking of the the phase space 
region |A| S [Vl — w^, 1] containing solitons, with the 
consequent compression of soliton eigenvalues on spectral 
gap edges at A = ±1. Such process leads to the collapse 
of the invariant measure towards two eigenvalues: 



limu(A)«-[<5(A 



1) + <5(A + 1)]. 



(52) 




FIG. 3. (Color Online). Time dependent dynamics of the 
many-body soliton ensemble for w = 1 [a), w = 0.95 (b), 
w = 0.9 {c), w = 0.85 (d). In all simulations the amplitude 
of the input -00 is set to p = 30. 



In this case the corresponding partition function scales 
as: 



lim 



(53) 



which predicts the emergence of a null free energy x = 0. 
The net effect of decreasing w is therefore to increase 
the " pressure" of soliton eigenvalues towards spectral gap 
edges, thereby leading to the collapse of the whole free 
energy landscape to a fiat profile. This picture well agrees 
with the numerical evaluation of the free energy x(/3,?«) 
(Fig. [2]). In particular, for w ^ Q the free energy shows a 
profile described by ()51|) . while for w = Q it dramatically 
collapses to zero. The geometric morphology of x(/3,u') 
can be then employed to derive a phase diagram for our 
thermodynamic system, with the shape of the free energy 
acting as an order parameter for different observable dy- 
namics. By considering Fig. [21 I identify two radically 
distinct dynamical evolutions: 

i) for (/3,w ^ 0) the free energy is of type ([5T|) . thereby 
supporting a first order phase transition in the soliton 
dynamics. The phase transition, born from the different 
scaling of the box-counting probability pi, leads to the 
accumulation of soliton particles with close but opposite 
eigenvalues (velocities) towards gap edges. Such a pro- 
cess, enforced by the action of the parameter w, is the 
signature of a strongly-discontinuous dynamics with soli- 
tons expected to immediately escape from their initial 
positions. 

ii) for (/3, w = 0), the free energy x collapses towards a 
flat profile x = 0- I^i this case no phase transitions are 
predicted, and no cooperative dynamics occurs. When 
ui = 0, in fact, the input profile -00 turns into the trivial 
constant solution of the NLS ?/'(x, t) = p, which does not 
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FIG. 4. (Color Online). Time position of the wave breaking 
point (dashed line) versus grayness w, with insets showing the 
corresponding field intensity profiles (solid lines). 



show any significant evolution. 

In order to verify this analysis, I investigate the time do- 
main dynamics of Eq. (1211) . Figures summarize the 
results for p = 30 and w varying in the range w € [0.9, 1]. 
In agreement with the phase diagram resulting from ([2|) , 
the system develops a discontinuous type of evolution 
known as dispersive shock [37','48'-5ll: the field initially 
collapses to a singular (wave breaking) point t* (Fig. [3] 
insets) and then generates a series of filaments escaping 
from t* (Fig. |4|). As discussed above, such dynamics 
stems from the many-body solitons cooperation result- 
ing from the first order phase transition supported by x- 
The phenomenon of eigenvalue compression, sustained by 
w, is then evident by comparing Figs. [S^-d. This effect 
can be then employed to control the soliton cooperative 
dynamics (Fig. [S]). As seen in Fig. [3l in fact, a change 
in w is accompanied by a significant shift of the shock 
point i*, which increases in time as soon as the grayness 
is reduced. The positive sign of the shift can be predicted 
from the free energy collapse at w — > 0. In this condition, 
^{x,t) = p and t* — cx), which anticipates an increasing 
of t* for a decreasing value of w. 

In summary, the thermodynamics analysis based on the 
soliton invariant measure predicts the emergence of a 
first-order phase transition, which gives rise to disper- 
sive shocks with completely controllable dynamics. Such 
a complex interaction, which cannot be captured with 
a reductionist analysis focused on the individual system 
component (i.e., the single soliton), is an emergent prop- 
erty of the whole soliton ensemble. 



IV. SOLITONS AND OTHER COMPLEX 
SYSTEMS 

As we have seen in the previous section, a many-body 
ensemble of solitons can act like a complex system, sup- 
porting the generation of emergent properties when co- 
operative dynamics among the soliton particles settles in. 



In this section, I will compare this system to various dy- 
namics manifesting complex phenomena. This analysis 
is not only important to highlight unique features of soli- 
tons, but also to identify common denominators at the 
origin of complex behaviors. In the following I will con- 
sider two main complexity classes, which also represent 
traditional examples in statistical mechanics: glassiness 
and anomalous scaling in time series dynamics. 
Glassiness is a complex phenomenon whose representa- 
tive examples, other than in real glasses, are found in 
a large number of contexts including error correcting 
codes , super cooled liquids, polymers and neural networks 
d, H, [52l. l53j] . Although the physics of the glass state is 
still under active research, different approaches based on 
mode coupling theory and/or replica-based mean field 
analysis highlight the same fundamental physics, even 
from different perspectives. Glassy systems, in particu- 
lar, show a high degree of frustration :54] in their ground 
state, which appears extremely degenerate and character- 
ized by an incredibly-large number of metastable states of 
equivalent energy. When the temperature is cooled down 
below a specific threshold, the complexity of the ground 
state manifests in the so called glass phase, characterized 
by an exponential number of energy minima separated 
by high potential barriers which localize the dynamics 
in the configuration space. Due to the huge number of 
energy traps, the system in the glass state explores only 
a subset of the available phase space. For this reason, 
glassiness could not be described within the ergodic hy- 
pothesis — which is at the basis of classical statistical 
mechanics — and needs a different approach. Breaking of 
classical thermodynamic assumptions is also observed in 
the statistical analysis of time series generated by com- 
plex dynamics, such as DNA sequences, earthquakes, so- 
lar flares and eye movements during spoken conversations 
(55l - l60| . In all these systems, time sequences are first con- 
verted into many distinct trajectories and then into prob- 
ability distributions. At equilibrium, the underlying sys- 
tem complexity yields probability functions with anoma- 
lous scalings, which cannot be interpreted with classical 
statistical mechanics arguments and whose origin is still 
largely debated ^ . 

In analogy with these systems, solitons complex dynam- 
ics are here observed due to the breaking of conventional 
thermodynamics rules. For the soliton systems, how- 
ever, breaking is not driven by frustration — such as in 
glasses — but conversely by the lack of any entropy con- 
straint in the soliton invariant measure. The applica- 
tion of a standard canonical measure on a purely kinetic 
system like solitons, in fact, is a well known results of 
classical thermodynamics and yields a trivial free energy 
profile with no cooperative dynamics jSj. The lack of a 
canonical measure, or equivalently, the absence of any en- 
tropy constraint in the soliton invariant measure, is then 
the key for the appearance of phase transitions and soli- 
ton cooperation effects. Such a "canonical breaking" is a 
peculiar property of solitons: in glassy systems, for exam- 
ple, the dynamics within a single energy trap is always 
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described by employing a reduced ergodic assumption, 
with canonical probabilities satisfying the entropy maxi- 
mization principle (see e.g., for more details). More- 
over, due to the complete factorization of the Hamilto- 
nian into single-soliton contributions [see Eq. no 
frustration exist in the soliton ensemble which exhibits a 
specific crystalline ground state with all eigenvalues con- 
centrated in a region of the phase as small as possible, 
where the single-soliton Hamiltonian gets its minimum. 
In summary, the breaking of some assumptions of clas- 
sical statistical mechanics unifies different complexity 
classes, including solitons, glasses and complex time se- 
ries. Among them, a distinctive features of solitons is to 
posses an invariant measure which does not result from 
any entropy maximization principle. 



theoretical and experimental work aimed at discovering 
new emergent properties in the broad field of complex 
many-body dynamics. 
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Appendix A: Calculation of the spectral coefficients 

Cn 



V. CONCLUSIONS 

I addressed the problem of finding a thermodynamic 
description of a microscopic mechanical system described 
by an invariant measure which is not constrained by 
any a-priori condition. By taking a soliton ensemble 
as a relevant, interesting case, I demonstrated how 
to construct an invariant measure from the spectral 
transform of the system, and how to develop a ther- 
modynamics by employing Lie transformation groups. 
I have discussed a specific example with a universal 
model of wave propagation, namely the NLS equation, 
with reference to a transparent potential supporting an 
ensemble of gray solitons. For this system I showed, 
through detail calculations verified by a time-domain 
analysis, the existence of emergent properties which can 
be completely controlled by adjusting solitons input 
parameters. By a further comparison with different 
complex systems, I have discussed similarities and pecu- 
liarities of soliton waves. In particular, like glasses and 
complex time series, solitons break some assumptions of 
classical statistical mechanics and this is at the basis of 
the manifestation of their emergent properties. In the 
case solitons, in particular, breaking occurs because of 
the lack of any entropy maximization condition in their 
invariant measure, and this is a peculiar property of this 
ensemble. This work is expected to stimulate further 



In the numerator of c„, for A±„, I expand 5(A±„) in 
terms of trigonometric functions: 

6(A,„)= .."7" , =^, n = 0,...,M. (Al) 
I sm7r(p — n) i 

I then use de I'Hopital rule: 



^ =lim^-^^- 



a{sj) s->sj a{s) 



(A2) 



and then expand the spectral coefficient a in terms of 
complex fractions: 



S.J - S±r, 



(A3) 



with s±n — X±n + 2i(p — n) and n — 0, . . . ,M. By 
substituting Eq. (|X3)) into I obtain: 



(A4) 



1 n±„(sj-s±n) 



Equation (jAip . together (IA4[) . yields the searched expres- 
sion of c„ : 
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